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The energy conditions play an important role in the description of some important properties of 
the Universe, including the current accelerating expansion phase and the possible recent phase of 
super-acceleration. In a recent work we have provided a detailed study of the energy conditions for 
the recent past by deriving bounds from energy conditions and by making the confrontation of the 
bounds with supernovae data. Here, we extend and update these results in two different ways. First, 
by carrying out a new statistical analysis for q(z) estimates needed for the confrontation between 
the bounds and supernovae data. Second, by providing a new picture of the energy conditions 
fulfillment and violation in the light of the recently compiled Union set of 307 type la supernovae 
and by using two different statistical approaches. 
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I. INTRODUCTION 

In the study of the classical energy conditions [l[ in 
cosmological context, an important viewpoint is the con- 
frontation of their predictions with the observational 
data. Since the pioneering papers by Visser || a num- 
ber of articles have been published concerning this con- 
frontation by using model-independent energy-conditions 
integrated bounds on the cosmological observables such 
as the distance modulus and lookback time - E3 (see 
also the related Refs. [Til] ) . Energy conditions constraints 
on modified gravity models, such as the so-called f(R)- 
gravity, have also been investigated in Ref. [12j | and more 
recently in Ref. fl3| . 

In a recent work [lcj, we have shown that the fulfill- 
ment (or the violation) of these integrated bounds at a 
given redshift z is not sufficient (nor necessary) to ensure 
the fulfillment (or the violation) of the energy conditions 
at z. This amount to saying that the local confrontation 
between the prediction of the integrated bounds and ob- 
servational data is not sufficient to draw conclusions on 
the fulfillment (or violation) of the energy conditions at z. 
This crucial drawback in the confrontation between inte- 
grated bounds and cosmological data has been overcome 
in Ref. [lol ] , where new non-integrated bounds have been 
derived, and confronted with type la supernovae (SNe 



la) data of the gold lj| and combined [la ] samples. 

In this letter, to proceed further with the investigation 
of the interrelation between energy conditions on scales 
relevant for cosmology and observational data, we extend 
and update the results of Ref. [l(| in two different ways. 
First, carry out a new statistical analysis for q(z) esti- 
mates necessary for the confrontation between the non- 
integrated bounds and supernovae data. Second, we give 
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a new picture of the energy conditions fulfillment and 
violation for recent past (z < 1) by using the recently 
compiled Union sample [16[ with 307 type la supernovae 
along with the new as well as the previous [10j statistical 
tools. 



II. NON-INTEGRATED BOUNDS FROM THE 
ENERGY CONDITIONS 

In order to use the energy conditions on cosmological 
scale, we consider the standard cosmological approach 
in which the Universe is modelled by a 4— dimensional 
space-time manifold endowed with a locally homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric 



ds 2 



dt 2 -a 2 (t) 



dr 2 



1 - kr 2 



+ r 2 (d9 2 + 



sin 



, (1) 



where the spatial curvature k — 0,1 or —1 and a(t) is 
the scale factor. We additionally assume that the large 
scale structure of the Universe is determined by the grav- 
itational interaction, and hence can be described by the 
General Relativity theory. These assumptions restrict 
the energy-momentum tensor to that of a perfect fluid of 
density p and pressure p, i.e., T^ v = (p + p) u ll u v -pg^ v . 
In this context, the energy conditions take the following 
forms: [l[ 
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where NEC, WEC, SEC and DEC correspond, respec- 
tively, to the null, weak, strong, and dominant energy 
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conditions, and the density p and pressure p of the cos- 
mological fluid are given by 



the redshifts and distance modulus 
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where overdots denote the derivative with respect to the 
time t and G is the Newton's gravitational constant. 

The non-integrated bounds from energy conditions fiol ] 
can then be obtained in terms of the deceleration param- 
eter g(z) = —d/aH 2 , the normalized Hubble function 
E(z) — H(z)/Hq, and the curvature density parameter 
Ofco = —k/(a HQ) 2 , simply by substituting Eqs. §3§ into 
Eqs. © to give 
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where z = (ao/a) — 1 is the redshift, -ff(z) = a/a , and the 
subscript stands for present-day quantities. Here and 
in what follows we have used the notation of Ref. [Iol |. 
in which NEC, WEC, SEC and DEC correspond, re- 
spectively, to p + p > 0, p > 0, p + 3p > and p — p> 0. 

We note that for a given spatial curvature flko, NEC 
[Eq.©] and DEC [Eq. ©] provide, respectively, a lower 
and an upper bound on the E(z)—q(z) plane for any fixed 
redshift z*. The WEC bound [Eq.©] only restricts the 
normalized Hubble function for a fixed value ofQ^Q, while 
the SEC bound [Eq.©] does not depend on the value of 
the spatial curvature. Thus, for any given value of Qko, 
having estimates of g(z*) and E(z+) for different redshifts 
z*, one can test the fulfillment or violation of the energy 
conditions at each z ± . 

In this work, we focus on the FLRW flat (flko = 0) 
universe, in which the NEC and DEC bounds reduce, 
respectively, to q(z) > —1 and q(z) < 2. Now, the 
q{z it ) and E(z±) estimates are obtained by using a SNe 
la dataset, through a model-independent approach which 
consists in approximating the deceleration parameter 
q(z) function in terms of the following linear piecewise 
continuous function (linear spline): 
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where the subscript I means that the quantity is taken at 
zi , Azi = (z — z{) , and the prime denotes the derivative 
with respect to z. The supernovae observations provide 



p(z) = 5 log 
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where Sfe(x) = sin(a;), x, sinh(ir) for k — 1, 0, —1, respec- 
tively. Then, by using the following relation between q(z) 
and E{z): 

E(z) = exp /* 1±^£) dz , (11) 

along with Eq. (JTDJ) , we fitted the parameters of the q(z), 
as given by ©, by using the SNe la redshift-distance 
modulus data from the so-called Union sample as com- 
piled by Kowalski et al [16]. 



III. RESULTS AND DISCUSSIONS 

We have used two different statistical approaches to 
confront the energy conditions bounds with observational 
data. In the first approach, which holds only for the flat 
case 1 , we have computed the g(z*) estimates by marginal- 
izing over E(z+) and the other parameters (aj's) of the 
g(z) function, with la — 3d confidence levels (C.L.). In 
the second procedure (which was used in Ref. [13]), we 
have calculated the la — 3a confidence regions on the 
E{z- k ) — q{z- k ) plane, and used the upper and lower limits 
of q(z) to have the la — 3a C.L. of q(z) for all z < 1 
(recent past). 

To obtain a global picture of the violation and fulfill- 
ment of the energy conditions in the recent past by using 
the first statistical method, we have calculated the g(z*) 
estimates at 200 equally spaced redshifts in the interval 
(0,1], and our result are depicted in Fig. [1^, which shows 
the NEC, SEC, and DEC bounds along with the best- 
fit values and the la, 2a and 3a limits of q(z) in the 
q(z) — z plane. We note that WEC bound is fulfilled 
identically (E 2 (z) > 0) in the flat case. 

When an observational confrontation is needed for a 
non-fiat FLRW case (flko + 1 0) the second statistical pro- 
cedure has to be employed, because in this case the NEC 
[Eq. (JSJ] and DEC [Eq. ©] bounds depend on the esti- 
mates of E(z+). In this way, q(z+) estimates cannot be 
obtained by marginalizing over E(z*), and one has to cal- 
culate the confidence regions on the plane -E(z^) — g(z^) 
(second approach). In this work, however, we have used 
this approach also for the flat case (which we have fo- 
cussed on) in order to make a comparison of the obser- 
vational SNe confrontations obtained by using both sta- 
tistical procedures, and also with the results of Ref. [l(| • 



1 Note that in this case the NEC, SEC, and DEC bounds are 
independent of E(z). 
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FIG. 1: The best-fit, the upper and lower la, 2a and 3a limits of q(z) estimates, marginalizing over E(z) [panel (a)], and from 
the confidence regions on the E(z) — q(z) plane [panel (b)], for 200 equally spaced redshifts. The NEC and SEC lower bounds, 
and also the DEC upper bound for the flat case are shown. This figure shows that the SEC is violated with la confidence 
level until z ~ 0.4 [panel (a)] and z ~ 0.36 [panel (b)]. It also shows that, for both statistical methods employed, the DEC 
and NEC is violated for high redshifts within 3a confidence level, while the NEC is violated for z < 0.075 [panel (a)] and 
2 < 0.085 [panel (b)]. 



Figure [Bb) contains the result of our analysis obtained 
by using the second procedure. Besides the NEC, SEC, 
and DEC bounds, it shows the best-fit values, the up- 
per and lower la — 3a limits of q(z) from the confidence 
regions on the E(z*) — q(z+) plane calculated for each of 
previously used equally spaced 200 redshifts in the inter- 
val (0,1]. 

The two panels in Fig.Q]show the violation of the SEC 
with more than 3a in the redshift intervals (0.05, ~ 0.2) 
[panel (a)] and (0.07, ~ 0.19) [panel (b)], where the high- 
est evidence of SEC violation is at z = 0.133 with 5.69a 
[panel (a)] and 5.34c [panel (b)]. Unlike the result from 
the confidence regions approach displayed in panel (b), 
which indicates the breakdown of SEC within la — 3a in 
the whole redshift interval, we note that from panel (a) 
one has that, within la, the SEC is fulfilled for z > 0.64 
and is violated for z < 0.4. This indicates that, within 
la C.L., the Universe crosses over from a decelerated to 
an accelerated expansion phase for a redshift within the 
interval (~ 0.4, ~ 0.64). 

Regarding the NEC, Fig. Q] indicates its breakdown 
within 3(7 for low redshift intervals, i.e., (0,0.075) [panel 
(a)] and (0,0.085) [panel (b)]. For higher values of red- 
shifit, NEC is violated with 3a for (z > 0.89) [see panel 

(a) ] and for z > 0.825 [cf. panel (b)]. 

Concerning the DEC, wc note that its fulfillment takes 
place in most of the redshift interval for both statistical 
analyses [see panel (a) and panel (b)], but it is violated 
within 3a for z > 0.765 [panel (a)] and z > 0.74 [panel 

(b) ], which are intervals where the error bars of our esti- 
mates grow significantly. 

The comparison of the results obtained through the 



confidence regions approach by using the Union set of 
307 SNe la with those of Ref. [10( calculated through the 
same statistical procedure but by employing the 182 and 
the 192 supernovae of the gold and combined samples, 
shows that the errors bars of q(z ir ) from the Union sample 
analysis are smaller for redshifts lying in (0, ~0.7). The 
Union sample results reinforce the indication of the SEC 
violation and NEC fulfillment at low redshift pointed out 
recently in Ref. [l(| ■ 

Here, similarly to the analyses of Ref. we have 

found that the results of the analyses for the best fit, the 
upper and lower la values of Qko = — 0.00461°'™^ as 
given by five-year WMAP (l7| . are essentially the same 
of the flat case, with differences much smaller than the 
associated errors. We note that for the upper la limit of 
n k0 = 0.002, the WEC bound [E 2 {z+) > 0.002(1 + z„) 2 ] 
is fulfilled with 3a confidence level in the redshift interval 
(0, 1], while for the Q k0 interval (-0.0113,0) the WEC 
is identically satisfied. 

IV. CONCLUDING REMARKS 

In a previous work [Io[ we provided a picture of the vi- 
olation and fulfillment of the energy conditions in the re- 
cent past by deriving non-integrated bounds from energy 
conditions in terms of the deceleration parameter and the 
normalized Hubble function in the context of FLRW cos- 
mology, and made the confrontation of the bounds with 
SNe la data through estimates of q(z) from la — 3a con- 
fidence regions on the plane E{z) — q{z) calculated with 
gold [lij and combined [l5T | samples. 
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Here, we have extended and updated the confrontation 
between the non-integrated bounds and supernovae data. 
First, by using the fact that, in the flat case, the NEC, 
SEC, and DEC bounds do not dependent on E(z), we 
have carried out a new statistical analysis in which the 
q(z) estimates are obtained by marginalizing over the pa- 
rameters E(z) along with the q'^s of q(z) [see Eq. JSJ)]. 
Second, we have updated the previous work [1 0f | by pro- 
viding a new picture [see Fig. QJa) and [Fig. D^b)] °f the 
energy conditions fulfillment and violation from the re- 
cently compiled Union set of 307 SNe la along with two 
different statistical tools. 

On general grounds, our analyses indicate a possible re- 
cent phase of super-acceleration in which the NEC is vi- 
olated within 3ct confidence level for z < 0.075 [Fig. Ufa)] 
and z < 0.085 [Fig. QJb)] , and that the DEC is fulfilled 
with 3tr in the redshift interval (0,0.765) [Fig. Ufa)] and 
(0,0.74) [Fig. Hp))]. Regarding the SEC, our analyses 
show that, for both statistical approaches employed, the 
best-fit curve of q(z) crosses the SEC-bound curve at 



z i2 0.51, and that SEC is violated with 3cr within small 
low redshift intervals [Fig. QJa) and [Fig. DJb)] . 

Finally, an interesting fact that comes out of our SEC 
analysis with la C.L., obtained by using the recent SNc 
la Union set, is that for the new q(z) estimate [calculated 
by marginalizing over E(z)\ the deceleration to accelera- 
tion transition expansion phase of the universe took place 
in the redshift interval (~ 0.4, ~ 0.64 ). 
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